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■ Summary. The iron mass in galaxy clusters is about 6 times larger than could 

have been produced by core-collapse SNe, assuming the stars in cluster galaxies 
formed with a standard IMF. Type-la SNe have been proposed as the alternative 
dominant iron source. We use our HST measurements of the cluster SN-Ia rate at 
high redshift to study the cluster iron enrichment scenario. The measurements can 
constrain the star-formation epoch and the SN-Ia progenitor models via the mean 
qq ' delay time between the formation of a stellar population and the explosion of some 

\Q , of its members as SNe-Ia. The low observed rate of cluster SNe-Ia at z ~ 1 pushes 

£ — ■ back the star-formation epoch in clusters to z > 2, and implies a short delay time. 

' We also show a related analysis for high-z field SNe which implies, under some 

conditions, a long SN-Ia delay time. Thus, cluster enrichment by core-collapse SNe 
from a top-heavy IMF may remain the only viable option. 

1 The SN la Rate in z < 1 Galaxy Clusters and the 
O ' Source of Cluster Iron 

The iron mass in galaxy clusters is about 6 times larger than could have been 
produced by core-collapse supernovae (SNe), assuming the stars in the cluster formed 
with a standard initial mass function (IMF; e.g., Renzini 1997). SNe la have been 
k>( i proposed as the alternative dominant iron source. Different SN la progenitor models 

. predict different "delay functions" , between the formation of a stellar population and 

the explosion of some of its members as SNe la. We use updated measurements of the 
total iron mass-to-light ratio in rich clusters to normalize the predicted SN la rate 
in clusters vs. redshift, using the delay function parameterization of Madau, Delia 
Valle, & Panagia (1998). We then use our previous measurements of the cluster 
SN la rate at high redshift (Gal- Yam, Maoz, & Sharon 2002) to constrain SN la 
progenitor models and the star-formation epoch in clusters. The low observed rate 
of cluster SNe la at z ~ — 1 (Fig. 1) means that, if SNe la produced the observed 
amount of iron, they must have exploded at even higher z. This puts a > 95% 
upper limit on the mean SN la delay time of r < 2 Gyr (< 5 Gyr) if the stars in 
clusters formed at Zf < 2 (z/ < 3), assuming H a — 70 km s _1 Mpc -1 (see Maoz & 
Gal-yam 2003 for full details). In the next section, we show that, for some current 
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versions of cosmic (field) star formation history (SFH), observations of field SNe la 
place a lower bound on the delay time, r > 3 Gyr. If these SFHs are confirmed, the 
entire range of r will be ruled out. Cluster enrichment by core-collapse SNe from a 
top-heavy IMF will then remain the only viable option. 
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Fig. 1. Predicted SN la rates vs redshift, if most of the iron mass in clusters is 
produced by type-la SNe following a brief burst of star formation at redshift Zf = 2 
(solid curves) and Zf — 3 (dashed curves). The bottom panel is a zoom on the 
low-z region of the top panel. The different curves are for SN la transfer functions 
with mean delay times, r, as marked. Cluster SN la rate measurements are by Reiss 
(2000) and Gal- Yam et al. (2002). The latter are shown with 95%-confidence vertical 
error bars. The horizontal error bars give the visibility-time-weighted redshift ranges 
of the cluster samples. The z; — 2 models (solid curves) with r > 2 Gyr are clearly 
ruled out by the z ~ 1 SN-rate measurement, even after accounting for a 30% 
uncertainty in the nomalization of the models. The zj = 3 model (dashed curves) 
with t = 5 Gyr predicts unacceptably high rates at low z. 
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2 The Redshift Distribution of Field SNe la: Constraints 
on Progenitors and Cosmic Star Formation History 

In this section, we use the redshift distribution of SNe la discovered by the Supernova 
Cosmology Project (Pain et al. 2002) to constrain the star formation history (SFH) 
of the Universe and SN la progenitor models. Fig. 2 illustrates how, for a given choice 
of parameters describing the SFH and SN la delay time, we predict the observed 
cumulative (i.e., unbinned) redshift distribution in a given survey, and compare it to 
the data. Given some of the recent determinations of the SFH, the observed SN la 
redshift distribution indicates a long (<; lh^ 1 Gyr) mean delay time between the 
formation of a stellar population and the explosion of some of its members as SNe la 
(Fig. 3). For example, if the Madau et al. (1998) SFH is assumed, the delay time 
r is constrained to be r > 1.7(r > 0.7)ft _1 Gyr at the 95%(99%) confidence level 
(CL). SFHs that rise at high redshift, similar to those advocated by Lanzetta et al. 
(2002), are inconsistent with the data at the 95% CL unless r > 2.5h~ 1 Gyr. Long 
time delays disfavor progenitor models such as edge-lit detonation of a white dwarf 
accreting from a giant donor, and the carbon core ignition of a white dwarf passing 
the Chandrasekhar mass due to accretion from a subgiant (e.g., Yungelson & Livio 
2000). The SN la delay may be shorter, thereby relaxing some of these constraints, 
if the field star formation rate falls, between z — 1 and the present, less sharply than 
implied, e.g., by the original Madau plot. The discovery of larger samples of high-z 
SNe la by forthcoming observational projects should yield strong constraints on the 
progenitor models and the SFH (see Gal- Yam & Maoz 2003 for full details). In the 
previous section we have demonstrated that if SNe la produce most of the iron in 
galaxy clusters, and the stars in clusters formed at z ~ 2, the SN la delay time must 
be lower than 2 Gyr. If so, then the Lanzetta et al. (2002) SFH is inconsistent with 
the data presented here (Figure 3) . 



3 Conclusions 

The large mass of iron in clusters, combined with the low SN la rate we have mea- 
sured in z ~ 1 clusters, require that, if SNe la are the dominant iron source, their 
delay time must be short. This would cast doubt, e.g., on the double degenerate 
models. The observed redshift distribution of field SNe la implies a slowly chang- 
ing field SNR(z) for SNe la. For several popular star-formation histories, this then 
indicates a long delay time. Taken together, these constraints may suggest that the 
iron in clusters is from core-collapse SNe, from an early stellar population with a 
top-heavy IMF. 
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Fig. 2. Illustration of the modeling and comparison to data. The SFH (&(z)) is 
modelled by two power laws, smoothly joined at z = 1.2. We denote the high-z index 
with q and the low-z index with (3. Panel a shows two examples, a "Madau" SFH, 
with a peak at z = 1.2 (solid curve), and a shallower model (dashed curve) reflecting 
the proposed modifications by Cowie et al. (1999) and Steidel et al. (1999). Panel 
b shows two examples of the expected SN la rate density following a brief burst 
of star formation. These delay functions are calculated using the prescription of 
Madau, Delia Valle, & Panagia (1998), with characteristic exponential delay times 
of t = lh^ 1 Gyr (solid) and r = 3/jT 1 Gyr (dashed). For display purposes, an 
arbitrary redshift of z = 3 has been chosen for the burst of star formation. SFH 
models are convolved with a delay function, and the resulting SN rate densities 
n(z) for a "Madau" SFH with r = 1ft" 1 Gyr (solid) and a "Cowie-Steidel" SFH 
with r = 3/i _1 Gyr (dashed) are shown in panel c. Panel d shows the predicted SN 
distributions, N(z), for the models of panel c, in a survey with the same observational 
parameters of the SCP search. KS tests show that the cumulative version of N(z) 
from a model combining a "Madau" SFH with a typical delay time of t = lh^ 1 
Gyr (panel e) is ruled out by the data, while a model with "Cowie-Steidel" SFH 
and t — 3/i _1 Gyr is consistent with the data (panel f). This "forward modelling" 
type of analysis is more powerful than the derivation of an observed SN rate at a 
mean z from the data (e.g., Pain et al. 2002; Tonry et al. 2003) and comparison to 
predictions (as in panel c), since no binning over redshift is carried out. Vertical axis 
units are arbitrary in panels a-d. 
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Fig. 3. Probability of SN la time delay values r, given the data, for particular SFH 
models. Assuming the SFH models shown in the upper panels, we can constrain the 
allowed values of r by the probability derived from the KS test (lower panels) . Points 
below the upper and lower dotted lines are ruled out at 95% and 99% confidence, 
respectively. 
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